Dynamical symmetry breaking in the 
Nambu-Jona-Lasinio model with external 
gravitational and constant electric fields 

E. Elizalde 1 ' 2 ': Yu. I. ShiPnov 1 ' 3 ' 1 " and V. V. Chitov 3 '* 

1 Consejo Superior de Investigaciones Cientificas, 
IEEC, Edifici Nexus-204, Gran Capita 2-4, 08034, Barcelona, Spain 
2 Department ECM, Faculty of Physics, University of Barcelona, 
Diagonal 647, 08028, Barcelona, Spain 

3 Department of Theoretical Physics, Faculty of Physics, 
Kharkov State University, Svobody Sq. 4, 310077, Kharkov, Ukraine 



Abstract 

An investigation of the Nambu-Jona-Lasino model with external constant elec- 
tric and weak gravitational fields is carried out in three- and four-dimensional 
spacetimes. The effective potential of the composite bifermionic fields is calculated 
keeping terms linear in the curvature, while the electric field effect is treated ex- 
actly by means of the proper-time formalism. A rich dynamical symmetry breaking 
pattern, accompanied by phase transitions which are ruled, independently, by both 
the curvature and the electric field strength is found. Numerical simulations of the 
transitions are presented. 
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1 Introduction. 



The Nambu-Jona-Lasinio (NJL) model jl] has long been considered one of the most 
suitable approximations providing a correct description of low-energy strong interactions 
physics. It has been seen that dynamical symmetry breaking (DSB)[] and dynamical 
fermions mass generation take place in the models with four-fermion interactions, those 
of NJL and Gross- Neveu ||. 

External conditions, such as a non-zero temperature, a finite chemical potential or 
classical gauge fields, have been shown to enrich the phase structure of the NJL model 



essentially (for reviews and complete lists of references see, for instance, pl-pOl) 



The influence of external magnetic and electric fields in these models has been studied 



for some time [0]-[2O|. It has been found out that a non-zero magnetic field always breaks 



the chiral symmetry, while the presence of an electric field (EF) tries to restore it. 

Investigations of the influence of a classical gravitational field on the DSB phenomenon 
in the NJL model have been carried out for some years. It has been shown that curvature- 
induced phase transitions exist and that they must be incorporated into any realistic sce- 



nario of the early Universe [pl]|-||24| (for a general introduction to quantum field theory in 
curved spacetime see [ESI and for a recent review of the NJL model in curved spacetime, 
PB[). It turns out that, in spite of the relatively small value of the curvature-dependent 
corrections at the low energy scale to be investigated within the NJL model, these correc- 
tions appear to be inescapable, in the sense that they must be taken into account when 



one performs the necessary "fine tuning" of the different cosmo logical parameters [27 



On the other hand, an external electromagnetic field plays an important role in the 
early Universe, which both can contain primodial magnetic fields f28| , |29[| and have a 
very large electrical conductivity [^S, 30|. Investigations of the influence of an external 
magnetic field within the DSB in curved spacetime have been carried out recently pl| - 
p3| . Therefore, it seems natural to study the effect of a constant EF upon the DSB 



phenomenon in a curved spacetime extending the results of fll2 , I3[ , which have been 
done for the case of flat spacetime. 

In this paper we investigate the behaviour of the NJL model in the presence of an 
external constant EF treated nonperturbatively in the proper-time formalism ||. The 
linear-curvature corrections to the effective potential (EP) of composite bifermionic fields 
are calculated and the phase structure of the model is investigated for negative values 
of the coupling constant. Both the four- and threedimensional cases are discussed. It is 
well known that in an external EF particle creation takes place and that the EP acquires 
an imaginary part @-[13[. That is why we study the small EF limit when the particle 
creation velocity is negligible. The phase transitions accompanying the DSB process on 



J This phenomenon and its application in high energy physics have been described in 
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the spacetime curvature, as well as the values of EF strength will be described numerically 
in detail. 



2 The fermion Green function in an external con- 
stant electromagnetic field with linear- curvature 
accuracy. 

Our starting point is the Nambu-Jona-Lasinio model in curved space-time of arbitrary 
dimension, d, as given by the following action M : 



S 



d d xyf^g \ itp'j fl (x)D fl '4j + 



2N 



where the covariant derivative Dn includes the electromagnetic potential A p : 



Df, = 8^ - ieA^ + -oo ab ll a ab . 



(2) 



The local Dirac matrices Jn(x) are expressed through the usual flat ones 7 a and the 



tetrads e" ;: 



The spin connection has the form 

1 



ab 



(Tab = ^[7a,7fe]- 



1 



(3) 
(4) 



e av {d„e b v - d u e b ) + -e ai ' e b ?e Cil {d p e c v - d u e 



\e hv {d,el - d u e«) - ± e b »e a »e cll (d p e c u - d u e 



i, i , ■ <■// . ■■-/,' ,, p i ■ (5) 

where N is the number of bispinor fields ip a . The spinor representation dimension is 
supposed to be four. Greek and Latin indices correspond to the curved and flat tangent 
spacetimes, respectively. 

Introducing the auxiliary fields: 



A 

'N 



we can rewrite the action (1) as: 



S 



N 



d d xy/=g \ iifrfDpip - —(a 2 + it 2 ) - tfj(a + iirj^ij) 



Then, the effective action in the large- N expansion is given by: 



—T eff (a,ir) = - J d d x^g ° ^ - zlndet {i^(x)D tl - (a + ijs^)} 



(6) 



(7) 



(8) 



Here we can put ir = 0, because the final expression will depend on the combination 
er 2 + it 2 only. This means that we are actually considering the Gross- Neveu model [[| . 

Defining the EP as V e ff = —T e ff/N J d d x^/—g, for constant configurations of a and 7r 
we obtain: 



V eff = — + *Sp ln(x\[i^(x)D„-a}\x) 
By means of the usual Green function (GF), which obeys the equation 

(27 M L> M - cr) x G(x, x', cr) = S(x- x'), 

we obtain the following formula 

V^ff(cr) = — — iSpG(x,x,a). 



(9) 



(10) 



To calculate the liner curvature corrections the local momentum expansion formalism is 
the most convenient one [53|. In the special Riemannian normal coordinate framework 
we have: 

(12) 



= 5% + l -R% aa y p y% 



(13) 



V H^ah = -it fiXy (Tab, 



(14) 

y = x — x' . (15) 
The vector potential of the external electromagnetic field is chosen to be of the form: 

1 



A^x) 



-F r u 
2 v u ' 



(16) 



where F^ u is the constant matrix of the electromagnetic field strength tensor. Substituting 
Eqs. (12)-(16) into Eq. (10), we arrive at the following equation for the GF: 



G{x,x',a) = 5{x-x') (17) 

Fulfilling the expansion on the different spacetime curvature monomials 

G = G + G 1 + ..., (18) 

where Go is the GF in flat spacetime, G\ ~ R, and so on, we obtain the iterative sequence 
of equations: 



% ft + e 4L(x) — cr 



Gq(x, x) = 8(x — x) 



(19) 



% ft + e — ° 



Gx(x,x',a) + 
4 



-l a R% aa y p y a (d,-ieA^x)) 



i « a id „ ,A Jjc 



Gq(x, x', a) = 0. 



(20) 



Here and below we can already forget about the difference between the two kinds of 
indices, Greek and Latin, because it only shows up beyond the linear curvature approxi- 
mation in Eq. (20). We introduce the factor: 



x') = exp ie / A fJl (x")dx'^ l 
L J x' ' 

which is the solution of the equation 

{d„ - ieAJ x <f>(x, x') = 
and assume that, just as for flat spacetime, the GF has the form 

G(x, x\ cr) = $(x, x')G(x - x', a). 



(21) 



(22) 



(23) 



Then, the evident dependence on A^x) disappears and we obtain the equation deter- 
mining the GF Gi(x — x', a): 



a^Giix - x', a) = --7°i2Va2/V^G (x - x', a) 

o 



l a a bc R bcaX y x G (x-x',a) 



(24) 



The GF for the case of an external electromagnetic field in flat spacetime is supposed 
to be known. Denoting by G 00 (x — x', a) the GF that satisfies the equation: 



we get: 



Finally, 



x 



{i ft — ct)Goo(x — x' , a) = 5(x — x'), 



dx"G 00 (x — x" , a)Gi(x" — x 1 , a) 



---fR^aix - x'Y(x - xrd*G (x - x', a) 
o 

-Wo hc Rbca\{x ~ x') x G (x - x', a). 



G \(x — x , a) = J dx"Goo(x — x", a) 



—^1 a R%«a{x" - x'Y{x" - xTdfGoix" - x', a) 



(25) 



(26) 



-\i a V bC R b cax(x" -X'f 



G n (x" — x'a). 



(27) 
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We actually need the coincidence limit x — > x' in order to calculate the EP (11). 
This provides us with the opportunity to simplify (27), especially for constant curvature 
spacetime, where 



R„ 



R 



(2? 



W d{d-l) 

Thus, our basic expression for the GF Gi in a space-time of arbitrary dimension d, is the 
following: 



iR 



dzG 00 {-z,a) 



2^dM (z,a) 



12d(d- 1) 

-2z 2 Yd^G (z, a) + 3(d - 1) ^G (z, a)}. (29) 
This expression can be substituted into Eq. (11) directly, because <&(x,x) = 1. 



3 Effective potential and phase transitions in the 
NJL model in curved spacetime. 



The GF for the case of flat spacetime is found to be in the proper-time representation 

oo 

G (z,a) = e~ ind/i 



ds __^ stJ 2 , % /"iiiv 

{Ai\s) d l 2 4s 



es 



rcothr-yT^F^ 



G 00 (-z,a) = e 



-i-nd/A 



dt 



(Ant) d / [ 



-exp 



-i(<rt + 



At' 



1 

2t 



where: 



C ' nv v^v R/j. R> 



1 - (eEs) coth(e£s) 



E 2 



(30) 
(31) 

(32) 



Let us now consider the 3D case. After a Wick rotation is — > s, it — > t, we have the 
following expression: 



Sp £1(0,0-) = / dtds 



3 -(i+s)cr 2 



X 



72vr 3 / 2 7 (t + s) 3 / 2 (l + /tcotr) 2 
2k(k + t) + (9r + 5/t) cot r + k(t — 3k) cot 2 r 



(33) 



where r = eEs, k = eEt. Performing the integration over a, we get the EP: 



<7 



V eff (a) = — + 



2A 4tt 3 / 2 J s 5 /2 
1/A 2 



ds 



t cot r 



6 



R 



dsdte-^* 2 



x 



144vr 3 / 2 J J (t + s) 5 /2(i + /tco t T ) 2 

l/A 2 1/A 2 

-2k(k + r) + (9r + 5k) cot r + k(t — 3k) cot 2 r 



(34) 



This expression coincides exactly, in the limit E — > 0, with previous results |31| : 

.2 A 3 



a 



2A + 6^1 (l -4 )eXp( 4 ) + 2 ^ erfc( X : 



_R_ 

4A2 



exp(- 



o , 
A 2 ' 



7T— erfc( 



where 



erfc(a;) 



A 



A' 



(35) 



(36) 



There are two ways of justifying the introduction of the A parameter in equations (34) 
and (35) for the EP. The first one is the standard renormalization procedure, by means 
of the UV cutoff method. Then, in the limit A — ► oo, after the renormalization of the 
coupling constant 

I i \ 

(37) 



1 



1 

A 



7T^ 



we have the well known expression for the renormalized EP in flat spacetime [13, 14 



with linear curvature corrections [31 



V e E fJ» 



a 2 a 3 Ra 
2X R + 3vr + 24vr' 



(38) 



Three-dimensional four-fermions models have been shown to be renormalizable in the 
leading large- N order ||. Furthemore, external electromagnetic and gravitational fields 
do not interfere with the renormalization procedure, because the only divergence that 
appears in the expression for the EP (34) has already been removed by the substitution 
(37). This is due to the obvious fact that the local feature of renormalizability cannot 
be spoiled by the weak curvature of global spacetime or by an external EF. Formally 
this can be proved by looking at the calculations of the leading terms of the integrand 
in the limit s — > 0,t — > 0, which are the only essential ones in order to determine the 
UV-divergences of the EP (34). Thus, the renormalization has been performed by the 
formula (37) completely. 

However, in general, four-fermion models may be considered as low energy effective 
theories being derived from a more complete version of a quantum field theory (QCD, 
for example). In this case the parameter A can be treated as a natural characteristic 
scale limiting the range where our low-energy approximation is valid, and then our model 
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will describe some phenomenological effects of elementary particle physics. Bearing these 
reasons in mind, we can otherwise maintain A finite and investigate the DSB phenomenon 
in our model with this fixed cut-off. 

The integrand contains the function cot(e.E's), that periodically goes to oo. It developes 
an infinite set of poles on the integration path and obligues us to take into account the 
contribution of the corresponding residua, given by the imaginary part of the EP. The 
presence of imaginary terms in the EP means that particle creation takes places and that 
our vacuum is actually unstable. The solution to this problem lies outside the limits of our 
present investigation. The simplest possibility seems to be to consider a comparatively 
small EF strength, with specific values that would provide an exponentially depressed 
particle creation rate. 

To start with we shall consider the case of a 3D flat spacetime: 



oo 



1/A' 



\2 



where the integrand function has its poles at the points 



Till 

Sn = jE' 71 = 1 ' 2 '- (40) 



The standard residue technique yields: 

2 1 °° A 

1/A 2 

-*^4^X> 1 e M-^° % (41) 

where P. v. means "principal value." Using the analytical continuation of the Hurwitz 
^-function, we find the following expression for the renormalized EP: 



V effM°)- 2 \ R 4tt r C( 2>2teE } \2ieE) 



(42) 



A numerical analysis of KeV/J^ R (a) for negative coupling constant gives the typical 
behaviour of a first-order phase transition, as shown in figure 1. The critical values 
are defined as usual: E& corresponds to the strength of EF for which a local nonzero 
minimum appears, E c , when the real part of EP is equal at zero and at the local minimum, 
and E C 2, when the zero extremum becomes a maximum. For all figures, an arbitrary 
dimensional parameter, /x, defining a typical scale in the model, is introduced in order 
to perform the plots in terms of dimensionless variables. A obviously does not appear 
anywhere because after renormalization (37) it must be set such that A — > oo. 
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The curvature-dependent part of the effective potential does not contain new UV- 
divergences but, unfortunately, A must be kept finite here, because the integrals cannot be 
calculated in the general case. Therefore, an analytical continuation cannot be perfomed. 
However, as mentioned above, the limit of weak electric field is the most suitable one here 
in order to investigate DSB correctly. The most natural way is to keep the flat part of 
EP in the same nonpertubative form (42) meanwhile in the curvature-dependent terms 
which are small already by themselves only the weak EF limit is taken into account. 
Thus, neglecting the next order exponentially depressed terms in curvature corrections, 
the renormalized EP turns out to be the following: 



(3D) 



eff,R 



(2ieE) 



3/2 



2A 



R 



4:71 



2C(- 



2' 2ieE 



/ a 2 \V 2 1 
\2ieEJ 



Ra iR(eE) 



1/6 



+ 



24vr ' 2tt 2 3 7 /3 



■cxp( 



(43) 



The shape of ReV^j^j^cr) is shown in figure 2, for the case of a fixed curvature and 
coupling constant, and in figure 3, for the case of fixed EF strength and coupling constant. 
The shapes are absolutely typical of first-order phase transition pictures. 
For a finite cut-off scale A, the effective potential is given by: 



ReK 



(3D) 



eff 



( \ ° 

((7) = 2A 



47r 3/2 



P.V. 



ds 

7572* 



T COt T 



1/A 2 



R 



-P.v. 



dsdte-( t+s > 2 



x 



1447T 3 / 2 J J (t + s) 5 /2(l + KcotT )2 

1/A 2 1/A 2 V ' V ' 

-2k(k + r) + (9r + 5k) cot r + k(t — 3k) cot 2 r 



(eE) 



3/2 OO 



nn 



4vr 2 



n 3 ^ 2 exp( ; ;(1 



71=1 

oo 



dy- 



eE 

exp (£(<p-y)) 



R ( eE 1/2 ™ nn 2 

+ 144 l vr j eE ° ' J. 2 (! + V 2 ) 2 (y + nn- iff/ 2 



X 



a 



where 



2y(y 2 + 4) (7m - y? + y)— + 3{nn - if + y)(3 - y 2 ) + y(8 - y 2 ) 



ip = arccos • 



(44) 



(45) 



(46) 



In contrast to the case of D — 3, the 4D four-fermion models are not renormalizable. 
Thus, we have to keep A finite, and do just the same type of calculations as for the 
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previous 3D situation: 



SpG 1 (0,<r) = 



iRa 
96V 



dtds 



-(t+s)a 2 



(t + s) 2 (l + KCot r) 2 



x —k(k + r) + 2(k + 3r) cot r + 2k(t - n)cot r 
We arrive at the EP with a finite cut-off scale: 

^ 2 f 



(47) 



2A 8tt 2 



P.u. 



e S,T rcotr 



1927T 2 



l/A 2 

dsdte'^ 3 ^ 2 



(t + s) 3 (l + /tcot r) 2 



x 



l/A 2 l/A 2 

— k(k + r) + 2(« + 3r) cot r + 2/t(r — k) cot 2 r 

(4D) _ 



(48) 



(e£) 2 ~ 



8tt 3 



1927T 



71=1 
OO 



-2 / ^ 2\ 

exp(--<7 ) 



x 



eE/A 2 

a 2 



_|_ y 2 ) 2 (y -)- 7T71 — </?) J 



y(y +4)(7rn- </? + y) — + 6(vm- </?) - 2y(y -5) 



(49) 



However this cut-off regularization with a finite A is not the only one that can be 
performed. Instead of it, it is also possible to send A — > oo and simultaneously rewrite 
our formula for the EP in 4 — 2e dimension, introducing therefore some kind of dynamical 
regularization which seems to be more useful for numerical simulations. Then, for weak 
EF the EP can be expressed in the following, more convenient form: 



a 
2A 



{eEf 
8tt 2 



r(-i + e) 



4c( - 1+e '-^ )+ 4 



ifc 2 , iR(eE) 2 / 3 . (T 2 . . 2 2/ , 



(50) 



Figure 4 presents the real part of EP in the (4 — 2e) -dimensional case in flat spacetime. 
It has been found by the numerical analysis that the final results are almost independent 
of the special e values. The phase transition is of second-order, because there exists a 
unique critical EF strength, defined as the value for which the zero extremum becomes a 
maximum and, simultaneuosly, a nonzero absolute minimum appears. This critical value 
separates the situations where the zero extremum is a maximum and a non-zero absolute 
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minimum of the potential exists, from the ones where the origin is the only and absolute 
minimum. 

Figures 5 and 6 show the shapes of the real part of EP in curved (4 — 2e)-dimensional 
spacetime for fixed EF strength and spacetime curvature on the abscissa axis. In both 
cases, a second-order phase transition occurs. 

It should be emphasized that in the presence of an imaginary part in the EP, the 



criterion for restoration of the chiral symmetry must be modified |13| , PL8| . In fact, in 
that case the order parameter has to be chosen as a(ipip), but not simply as a because, 
although the later vanishes at the critical point, the former might still break the chiral 
symmetry through the corresponding term in the action (7). Formula (11) gives that 

N 

X 



(W) = iNlmVl ff (<j)--<T (51) 



and, therefore, we should only check whether our EP obeys the criterion: 

]imalmV' ff (a) = 0. (52) 

Formulae (43) and (50) make us sure that this criterion is satisfied here for both the 3D 
and 4D cases. 



4 Conclusions. 

Phenomenological models with four-fermion interaction seem to be useful for the de- 
scription of low- and inermediate-energy physics of strong interactions 0. The fermions 
are usually treated as quarks and the composite particles generated by the dynamical 
symmetry breaking as mesons, despite of some simplifications making the models under 
consideration more calculable. 

The nonzero vacuum expectation value of composite bifermionic operator (ipip) plays 
in DSB scheme just the same role as (if) in the traditional Higgs one and defines the 
dynamically generated mass of fermions in our case. 

Usually the symmetry to be broken under the DSB mechanism is the chiral one. A 
dynamical version of fermions mass generation and dynamical chiral symmetry breaking 
have been investigated very carefully and some fruitful applications for the real high- 
energy physics have been found. 

In particular, the NJL model in external electromagnetical field has been studied from 
different points of view to realize the crucial role of that field in the DSB for flat spacetime 
||- [2~0|. The effect of dynamical chiral symmetry breaking by external magnetic field 



at any attractive interaction of the fermions has been shown to be universal and model 
independent phenomenon | lq| . 



2 More detailed information and complete literature list can be found in |Q, for example 
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In our paper we have extended these investigations to the case of the presence of an 
electric field and non-zero spacetime curvature. Both these values try to restore the chiral 
symmetry in contrast to the magnetic field and their competition with the DSB catalysis 
by the magnetic field has to be understood more precisely. 

The Gross-Neveu model with the simplest four-fermionic interaction (ipip) 2 is invariant 
under the discrete chiral transformation 

if; -> 7 5 ^ (54) 

Because of the discrete type of this symmetry Goldstone modes do not appear. 
However NJL model has a continuous chiral invariance under the transformations 

$ _> e ^ (55) 

In contrast to the previous case Goldstone bosons have to appear if only the symmetry 
is being broken. In the phenomenological applications of NJL model 7r field is treated 
as Goldstone modes corresponding to the massless pions, which obtain the kinetic terms 
through the quantum corrections. Of course, if we take into account the quarks' current 
masses these modes become pseudo- Goldstone ones ]5|, We only split them out of our 
consideration from the very beginning by the appropriate choice of the direction in the 
a — it space and use this choice's advantage that DSB is ruled by a itself. 

The investigation of the phenomenon realization in the presence of magnetic field in 
flat spacetime has been performed both for NJL model and for quantum electrodynamics 
and quantum excitations of a and tt fields dispersion laws have been found [[HJ. To 
generalize this calculations for the curved spacetime case is an interesting problem we 
are dealing with now. 

However, here we have studied the phase structure of the NJL model in curved space- 
time with external constant pure electric field only. In three dimensions, a first-order 
phase transition takes place, both on the electric field strength and on the spacetime 
curvature. In the four- dimensional case - or, more precisely, in 4 — 2e dimensions with 
negligible e - the phase transition is a typical second-order one, on both of these external 
parameters. 

Our approximation is quite consistent, because all of the critical values are very small 
in comparison with the characteristic scale of the model fi. It should be noted that 
the consideration, in a nonperturbative way, of an external electric field provides very 
interesting effects even in flat spacetime, and that the linear curvature terms induce some 
additional processes of chiral symmetry restoration. 

We clearly observe that a positive spacetime curvature tries to restore chiral symmetry, 
as an external electric field does. This is in contrast with what happens in the case of 
a magnetic filed, where chiral symmetry is broken for any value of the field strength 
although positive curvature might restore chiral symmetry again for some critical value. 
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Figure 1: The behaviour in 3D of the real part of the effective potential Re V'^ R / ^ as a 
function of a/^i is depicted for R = and fixed X/j, = —100. The curves from the upper 
part to the bottom of the plot correspond to the following electric field strengths: eE/[i 2 = 
0.0003; 0.0023; 0.00020; 0.00017; 0.00005, respectively. The critical values, defined as usual, are 
given by: eE cl //j, 2 = 0.00023; eE c / V = 0.00017; eE c2 /^ 2 = 0.00005. 
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Figure 2: Behaviour in 3D of ReV e ff t R/fj, 3 as a function of a//i for fixed eE/ = 0.00005 and 
A/U = —100. From above to below, the curves in the plot correspond to the following values 
of R//J- 2 = 0.006; 0.005; 0.004; 0.0032; 0, respectively. The critical values, defined as usual, are 
given by: R cl /fi 2 = 0.005; R c /^ 2 = 0.0032; R c2 /^ 2 = 0. 
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Figure 3: Behaviour in 3D of Re Veff^/fi 3 as a function of cr//x, for fixed R/fi 2 = 0.0032 and 
A/i = —100. Starting from above, the curves correspond to the following values of eE/[i 2 = 
0.00015; 0.00011; 0.00007; 0.00005; 0.00001, respectively. The critical values, defined as usual, 
are given by: eE&/ '/i 2 = 0.00011; eE c /fi 2 = 0.00005, while eE&f/J, 2 does not exist for the given 
values of fi\ and R/ y?. 
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Figure 4: Behaviour of Re V e ff / fi^ 4 ~ 2t ^ in a (4 — 2e)-dimensional spacetime as a function of 
cr/fi, for fixed R = 0, e = 0.02 and A/U = —1000. Starting from above, the curves in the plot 
correspond to the following values of eE/fi 2 = 0.02; 0.01; 0.005; 0.001, respectively. The critical 
value is given by: eE c /fi 2 = 0.02. 
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Figure 5: Behaviour of Re V e ff / fi^ 4 ~ 2t ^ in a (4 — 2e)-dimensional spacetime as a function of 
cr/fi, for eE/fi 2 = 0.001, e = 0.005 and fixed \[i = —1000. Starting from above, the curves 
correspond to the following values of the curvature: R/fj, 2 = 0.0045; 0.0035; 0.02; 0.01; 0; -0.001, 
respectively. The critical value is obtained for Rc/fi 2 = 0.0035. 
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Figure 6: Behaviour of Re V e ff / fj,^ 4 ~ 2e ^ in a (4 — 2e)-dimensional spacetime, as a function 
of a/mu, depicted for R/[i 2 = 0.002, e = 0.005 and fixed = —1000. From above, the 
curves in the plot correspond to the following values of the electric field strength: eE/fi 2 = 
0.004; 0.0028; 0.02; 0.015; 0.001, respectively. The critical value is reached at eE c /^ 2 = 0.0028. 
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